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ABSTRACT 

With the use of iV-body calculations the amount and properties of escaping stars from 
low-iV (N = 100 and 1000) young embedded star clusters prior to gas expulsion are 
studied over the first 5 Myr of their existence. Besides the number of stars also different 
initial radii and binary populations are examined as well as virialised and collapsing 
clusters. It is found that these clusters can loose substantial amounts (up to 20%) of 
stars within 5 Myr with considerable velocities up to more than 100 km/s. Even with 
their mean velocities between 2 and 8 km/s these stars will still be travelling between 
2 and 30 pc during the 5 Myr. Therefore can large amounts of distributed stars in 
star-forming regions not necessarily be counted as evidence for the isolated formation 
of stars. 

Key words: stellar dynamics - methods: iV-body simulations - binaries: general - 
stars: formation - open clusters and associations: general 



1 INTRODUCTION 

Stars do not form in absolute isolation but rather in 
groups or embedded clus t ers in dense molecular cloud 
cores dLada fe Ladal 1 19951 . 120031 ; lAdams fc Mversl l200ll : 
I Allen et alj|2007t h Theoretical and observational results in- 
dicate that these groups and the majority of the em- 
bedded clusters dissolve quickly due to gas expulsion 
after about 10 Myr and release their stars into the 
galactic field jKroupa et al.ll200ll: iKroupa fc Bouvierll2003l; 



Adams fc Mversl l200ll; iLada, fc Ladal 120031 ; IWeidner et ail 
20071 ; iBaumgardt fc KroupalbOOTl h But the ratio of the dis- 
tributed mode of star formation and the clustered one is 
not fully understood. Observations show a distributed frac- 
tion of stars in th e Orion molecular clouds of about 20% 
jAllen et al.ll2007t ). It is generally assumed that these frac- 
tion is primordial as dynamical evolution of the clusters is 
considered to slow to reproduce the observed fraction of iso- 
lated stars at such a young stage. But for a full physical 
understanding of the formation of stars it is vital to know 
all modes of star formation in as much detail as possible. The 
purpose of this work is therefore to access whether or not 
the observed amount of distributed stars outside embedded 
clusters could have a dynamical origin and henceforth being 
initially formed in the young, tight embedded clusters or if 
there is intrinsically isolated star formation. Especially in 
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star-forming regions which do not show massive star clus- 
ters the amount of stars released through the dynamical 
evolution of low- AT clusters is not well known. In order to 
study this question a large series of A^-body calculations is 
performed and the number, the velocity and the mass spec- 
trum of stars ejected due to dynamical interactions within 
5 Myrs from this clusters are studied. 



2 THE CLUSTER SETUP 



With the use of the Af-body6 code (|Aarsethlll99 9a b, 2003, 
2008) the evolution of numerical star clusters with small 
numbers of stars are examined. Two different numbers of 
stars, N, are used, 100 and 1000. As the actual number of 
escaped stars can be rather low each set of initial conditions 
is simulated 100 times using different random number seeds 
to set up the masses, velocities and positions of the stars in 
order to improve the statistical significance of the results. 
All s tars in the clusters follow the canonical IMF (|Kroupal 
12002ft . 

Besides the initial number of stars several other initial 
parameters are studied with each 100 calculations. 

• All stars are initially single stars, 

• all stars reside in binaries, 

• 50% of the stars are in binaries, 

• the cluster is initially in virial equilibrium , 

• the cluster is initially sub- virial (collapsing with the ki- 
netical energy being one-tenth of the potential energy). 

• All star clusters are setup as Plummer-spheres 
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i|Plummerlll91ll ) with different initial half- mass radii of 0.1, 
0.25 or 0.5 pc. 

Each of these setups are then numerically evolved 
with A-body6 for 5 MyiQ. The evolution time of 5 Myr 
is chosen to minimise mass loss due to stellar evolu- 
tion, to avoid supernovae and primarily to have a set- 
ting of a star cluster still embedded in its parental 
cloud, prior to gas expulsion. Though, c lusters initially 
probably form from clumpy s tructures ( W illiams et al.l 
|2000| ; ICarpenter fc H odapp 2 0081). Simulations of sub-virial 
clumpy initial conditions ( Allison et al.l [2009h have shown 
that such clusters mass segregate faster than virial clusters 
due to the formation of short-lived very dense cores. But 
these calculations did not include a background gas poten- 
tial. 

To emulate the gas in which the stars are embedded dur- 
ing this time, t he cluster is setup with an additional plum - 
mer potential l)Plummerl Il91ll : IBinnev fc Tremainel Il987l ) 
of the same half-mass radius and the same mass as the 
cluster. Thus assuming a star-formation efficiency (SFE) 
of 50 %. Though, global SFE's are observationally around 
30% (|Lada fc Ladall2003l ), locally within the region where 
the cluster actually fo rms the value is most likely higher 
ijMoeckel fc Batdl2010h . A SFE of 50% has therefore been 
chosen. A large differences in results between 30 and 50% 
SFE is not to be expected. The gas potential is in all 
cases kept constant during the whole calculation time. 
As the results are qualitatively similar to A-body studies 
without a gas backgrou nd potential (e.g. iKroupal Il998l ; 
iKroupa fc Bouvierl 120031 ) , the exact shape of the potential 
does not seems to be overly important. But detailed further 
studies are needed to clarify the full impact of the shape and 
depth of the potential on the properties of escaping stars. 

The decision to include primordial binaries into the 
study is based on the observational r esults that open 
clusters host large number of binaries (|Sana et all 120091 ; 
ISollima et a l. 2010). The initial setup of the binary prop- 
erties increase the total parameter space of the calculations 
significantly. As a cluster of the richness of the Orion Nebula 
cluster stellar dynamics may change the mass function and 
bin ary properties in the core significantly in less than 10 6 
vr (iGoodwin fc Bastian 2006; Pflamm-Altcnburg fc Kroupal 
I2OO6I ; iPfalzner fc Olczakl 120071 : 1 Allison et al J 120091 ) . Andas 
the majority of the present-day field binaries are being pro- 
cessed through star clusters, their pro perties are not suited 
to be used for the nume rical setup (iDuauennov fc Mayor! 
Il99ll : lKroupdfl99^ ldlrblkl: IKroupa fc Bouvierll2003l ). There- 
fore, the results of the inverse dynamical population synthe- 
sis of iKroupal (|l995bl ) are used, which predict a flat (ther- 
mal) period distribution (resulting in a uniform logarith- 
mic semi-major axis distribution) and random p airing of 
the st ars in systems of low-mass (< 1 Mq) stars (|Kroupal 
Il995bl l3). While t here are indications for non-random pairing 
of massive stars l|Kobulnickv fc Frverj [20071 : IWeidner et al.l 
2009), the low-mass results are also used for massive bina- 
ries in this study. Firstly, because the mass where a change 
of properties might occur is not determined yet and it is 



1 Depending on the initial radius of the cluster 5 Myr are between 
about 10 and 150 crossing times. See eq. [2] for how to calculate 
the crossing time. 



not clear if mass-ratios near to unity for massive stars 
are primordial or alread y a sign of dynamical evolution 
(Pfalz ner fc Olc zak 2007) and secondly, because due to the 
choice of N = 100 and 1000, only a relatively small num- 
ber of massive stars is to be expected anyway. For example, 
there should be only three stars above 8 Mq for N = 1000 
and a Kroupa-IMF. 

It should be noted here, that rapid changes in the bi- 
nary properties might also at least partly due to the initial 
conditions of A-body calculations and do not really reflect 
the star formation process. I.e. that unstable binaries that 
are easily destroyed in early cluster dynamics may not even 
form in the first place l|Moeckel fc Batdl201ol ). 

In the N = 1000 cases, for the IMF a maximal mass, 
fttmax, of 25 Mq is chosen, in accord ance with the m max — 
M oc i-relation by I Weidner et all (|20ld ). In the N = 100 case 
this relation yields a m max of 7 Mq. But also comparison 
calculations are done using m max = 25 Mq for the N — 100 
case. 



3 RESULTS 

In the series of Figs. [1] to [8] the results of the A-body6 cal- 
culations after 5 Myr are presented. There the mean num- 
ber of lost stars, <Ai ost >, the mean lost fraction of mass, 
<Mi ost >/<Mtot>, the mean escape velocity, <v CBC >, the 
mean escape velocity dispersion, <a>, and the mean mass 
of the lost stars, <mi ost >, are studied in dependence of the 
initial cluster relaxation time, ireiax- Th e latter is calculated 
by us ing the following equations from IBinnev fc Tremainel 
l|l987l) : 



N 



8 In A' 



(1) 



with N being the number of stars and t CIOSS the crossing time 
a star needs to travel through the cluster which is given by, 
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with G being Newton's gravitational constant, _R oc i the clus- 
ter radius and M to t its mass. 

It is important to note that the relaxation time of a clus- 
ter is not a constant but also evolves with time. Though, for 
the clusters which start in virial equilibrium, the change is 
negligible over the short period of time studied here. The 
clusters expand about ~10% due to evolutionary and dy- 
namical mass loss. This leads to a ~10% increase in £ r eiax- 
But as the sub-virial clusters contract significantly during 
the calculation time span, their i rc iax decreases and becomes 
comparable to a cluster with half their half-mass radius. 

In the B panels of the Figs. [3] and [7] are addition- 
ally shown several literature descripti ons of star loss. The 
first o ne is the analytical description of IBinnev fc Tremainel 
l| 19871 . Eq. 8-85) based on the Fokker-Planck approximation 
for the mass loss by evaporation of stars from a cluster (with 
large A). 



M(t) = M 1 



7k c t 

relax 
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with Mo being the initial cluster mass, k c ~ 0.003 and ^ c i ax 
the initial relaxation time (eq. [I}. With the use of the mean 
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mass, m mean of the IMF, which is for a Kroupa IMF m mea n 
~ 0.36 Mq, the number of stars lost over time can be ap- 
proximated: 



N(t) = M(t)/m n 



(4) 



This equation is indicated with solid black dots in the two 
Figs. [3] and 

This description is strictly valid only for large N and 
does not seem to pr ovide a good characterisation of the 
iV-body calculations. iHeggid (|l974h assumed that in low-JV 
clusters binary interactions should dominate over the evapo- 
ration of stars. He formulated the following dissolution time 
scale for low-N clusters: 

- iV 2 

^diss — 100 ^ cross ' ^ ' 

which can also be used to estimate the star loss from a clus- 
ter in the following way: 



N diss = N e 



ln(JVo) 
*di B5 



(6) 



were No is the initial number of stars. This description is 
included as open squares in the Figs. [3] and [7] 

An third description of stellar loss can be derived from 
the evaporation time scale w hen assuming the collision less 
Boltzmann equation given in lBinnev fc Tremaind (|l987l ). 



= 136t rc la 



(7) 



which can be transformed into the following number loss 
formula: 



N e - 



N e 



ln("0) , 



(8) 



where iVo is the initial number of stars. In the Figs. [3] and [7] 
this relation is plotted as open triangles. 

The Figs. [1] to |4] show the results using the star clusters 
with initially N = 1000 after 5 Myr while the Figs. [S] to \E\ 
show the same for N = 100. 



4 DISCUSSION 

As is visible in Fig. [1] the number and mass fraction of lost 
stars for the N=1000 clusters is only dependent on the initial 
radius of the cluster (— > t re iax), while the binary fraction does 
not play a major Likewise in the case of a collapsing 

cluster (sub-virial, dashed-dotted line). A sub-virial cluster 
just behaves like a somewhat more concentrated cluster as 
they collapse and therefore change their timescales. 

In the case of the TV = 100 clusters (Fig. [5} there seems 
to be a dependence of the number loss on the binary fraction 
for objects with large radii (i rc iax > 2 Myr). But this might 
still be a low-number stochastic effect as only a couple of 
stars leave these clusters (N\ OBt < 5). When comparing the 
lost mass fraction of the N = 1000 and N = 100 calculations, 
the richer cluster seem to loose 2 or 3% more mass, reflecting 
the higher encounter rate in the richer objects. Allowing for 
more massive stars in the low-iV clusters (short-dashed lines 
in Figs. f5]to[7| increases the star loss by similar amounts as 
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Figure 1. Mean number of stars lost vs t re iax (Panel A) and mean 
percentage of cluster mass lost vs t re iax (Panel B). The solid line 
shows the case without any initial binaries, the dotted line is the 
50% initial binary case, the dashed line is the 100% initial binaries 
case and the dashed-dotted line marks the sub-virial case. 



2 It should be noted here that all figures are evaluated after a 
calculation time of 5 Myr and that the abscissae in most cases is 
the relaxation time and therefore not a time evolution. 
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Figure 2. Mean escape velocity, <t> C sc>, vs t re iax (Panel A) and 
mean velocity dispersion, <<r>, vs t re lax {Panel B). The solid line 
shows the case without any initial binaries, the dotted line is the 
50% initial binary case, the dashed line is the 100% initial binaries 
case and the dashed-dotted line marks the sub-virial case. 
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Number of relaxations during 5 Myr 

Figure 3. Panel A: Mean mass of the escaped stars vs t re i a x- 
The solid line shows the case without any initial binaries, the 
dotted line is the 50% initial binary case, the dashed line is the 
100% initial binaries case and the dashed-dotted line marks the 
sub-virial case. Panel B: Total number of lost stars vs number of 
relaxation times the cluster experienced during the 5 Myr sim- 
ulation time. The solid line shows the case without any initial 
binaries, the dotted line is the 50% initial binary case, the dashed 
line is the 100% initial binaries case and the dashed-dotted line 
marks the sub-virial case. The solid black dots are the predicted 
numbers from eq. [4] the open squares from eq. \E\ and the open 
triangles from eq. [8] 
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Figure 4. Median and maximal flight ranges from the cluster 
centre of the escaped stars vs t re i ax . The solid line shows the 
case without any initial binaries, the dotted line is the 50% initial 
binary case, the dashed line is the 100% initial binaries case and 
the dashed-dotted line marks the sub-virial case. 



a high-binary fraction. This is probably due to the higher 
average mass of the systems compared to single stars. A 
higher upper mass limit of the IMF likewise increases the 
average stellar mass. A larger average mass leads to higher 
momenta and therefore to larger velocities after 2- and 3- 
body interactions which are then more likely to exceed the 
escape velocity of the cluster. 

Like the mass loss, the mean escape velocity, <u OS c> 
(Panel A of Fig. [21, is mainly depended on t rc i ax , though the 
existence of a high-binary fraction raises <v csc > by about 
1 km/s, independent of £ rc iax- Primordial binaries introduce 
3- and 4-body interactions early-on in the cluster evolution 
opposed to the mainly 2-body interactions for single star- 
only clusters. Such higher order interactions generally result 
in larger velocities. 

The mean velocity of the ejected stars in the case of N 
— 1000 is 4 to 8 km/s (Fig. [2]) and the escaped stars travel 
on the average up to 10 to 30 pc (Fig. [3| within 5 Myr de- 
pended on the original size of the clustefl But individual 
stars can be up to 1 kpc away, as the highest escape velocity 
encountered for the N = 1000 clusters was about 260 km/s, 
though the typical range of escape velocities is between 2 
and 40 km/s. These v alues are in the same range as the re- 
sults o f lKroupall|l99cf ) who studied clusters with N = 400 in 
a Galactic tidal field but without a gas background poten- 
tial. A comparison with observations seem to c ontradict the 
relativ ely large number of lost stars found here. iTobin et al.l 
(2009) study the surroundings of the Orion Nebula cluster 
(ONC) and find that the majority of stars with velocities 
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Figure 5. Mean number of stars lost vs t re iax (Panel A) and 
mean percentage of cluster mass lost vs t rc i ax (Panel B). The 
solid line shows the case without any initial binaries, the dotted 
line is the 50% initial binary case, the long-dashed line is the 100% 
initial binaries case, the dashed-dotted line marks the sub-virial 
case and short-dashed line the case with m max = 25 Mq. 



3 A velocity of 1 km/s is roughly equal to 1 pc per Myr. 
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Figure 6. mean escape velocity vs t re l a x {Panel A) and mean 
velocity dispersion vs t re i ax (Panel B). The solid line shows the 
case without any initial binaries, the dotted line is the 50% initial 
binary case, the long-dashed line is the 100% initial binaries case, 
the dashed- dotted line marks the sub-virial case and short- dashed 
line the case with m max = 25 Mq . 
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Figure 7. Panel A: Mean mass of the escaped stars vs t re i ax - 
The solid line shows the case without any initial binaries, the 
dotted line is the 50% initial binary case, the long-dashed line is 
the 100% initial binaries case, the dashed-dotted line marks the 
sub-virial case and short-dashed line the case with m max = 25 
Mq. Panel B: Total number of lost stars vs number of relaxation 
times the cluster experienced during the 5 Myr simulation time. 
The solid line shows the case without any initial binaries, the 
dotted line is the 50% initial binary case, the long-dashed line is 
the 100% initial binaries case, the dashed-dotted line marks the 
sub-virial case and short-dashed line the case with m max = 25 
Mq. The solid black dots are the predicted numbers from eq. [4] 
the open squares from eq. [6] and the open triangles from eq. [8] 
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Figure 8. Median and maximal flight ranges from the cluster 
centre of the escaped stars vs t re i ax . The solid line shows the 
case without any initial binaries, the dotted line is the 50% initial 
binary case, the dashed line is the 100% initial binaries case, the 
dashed- dotted line marks the sub-virial case and short- dashed line 
the case with m ma x = 25 Mq. 



differences to the cluster larger than 10 km/s have little or 
no IR-excess and are therefore most likely old and not asso- 
ciated with the ONC. When looking only at escapees with 
velocities larger than 10 km/s after only 1 Myr (the age of 
the ONC. iHillenbrand fc Hartmann|[l998l ). only between 
and 14 are to be expected from the TV-body calc ulations, de- 
pendi ng on the initial conditions. Interestingly, iTobin et al.l 
(2009) do find 6 stars with a velocity offset above 10 km/s 
which have an IR-excess consistent of being young stars and 
therefore could have been ejecte d from the ONC, consis- 
tent with our predictions. As the ITobin et al.l (|2009l ) study 
can not differentiate between bound and unbound cluster 
stars below 10 km/s it is currently not possible to test the 
predictions of the TV-body calculations further. Also consis- 
tent with this picture is the recently discovered halo of low- 
ma ss stars around the small young cluster rj Chamaeleon- 
tis l|Murphv et al.|[2010h . Generally, though, is an IR-excess 
probably not a good age indicator for ejected stars as most 
of these had a close encounter with another star which lead 
to the ejection. During this encounter any disc around that 
star could be truncated or largely stripped, therefore reduc- 
ing or eliminating the IR-excess. In order to search for a 
distributed population of ejected young stars around clus- 
ters the position of the stars in the HR diagram will likely 
be a more robust age indicator. 

For less rich clusters (Panel A of Fig.|5J), smaller <v asc > 
are generally achieved. These are also increased by the pres- 
ence of binaries by 0.5 to 1 km/s. Stars escaping from TV — 
100 clusters on the average travel 3 to 6 pc (Fig. |SJ within 
5 Myr with a typical range of escape velocities between 1 
and 12 km/s. For the low-TV clusters the fastest escaper had 



about 40 km/s and travelled about 100 pc in total. Interest- 
ingly, the mean of the velocity dispersion of escape veloci- 
ties, a, is more strongly affected by the presence of binaries, 
both in the TV = 1000 (Panel B of Fig. [2} and TV = 100 
case (Panel B of Fig. [6]). Therefore, the presence of large 
fractions of binaries, while only slightly enlarging the mean 
escape velocity, adds more scatter to the velocities of the 
escaping stars. This increases the area around a young em- 
bedded cluster where to expect escaping stars significantly 
by 50 to 200%. The mean mass of the escaping stars, how- 
ever, seems to be rather independent of binaries in the TV 
= 1000 clusters but is strongly correlated with i rc i ax (Panel 
A of Fig. |3}. It rises from about 0.22 Mq in the 10 Myr 
case (— > -Red = 0.5 pc) to about 0.4 Mq for very compact 
(1 Myr — ¥ R cc \ = 0.1 pc) clusters, while for the collapsing 
clusters it stays at about 0.3 Mq. On the contrary for the 
low- TV = 100 clusters, <mi os t> roughly stays constant with 
^rciax (Panel A of Fig. [7} but varies strongly with the binary 
fraction. Here the sub-virial, 50% binary and higher upper 
mass limit case are rather similar and well separated from 
the only-single star case and the 100% binaries. 

In appendix|X]histograms of the mass and <v csc > of the 
escaped stars for the different initial conditions are shown. 
Mass histograms are only plotted for the cases without ini- 
tial binaries for TV = 1000 and TV = 100 (Panels A of Figs. [ATI 
and IA2|) . The mass histograms in the TV = 1000 look very 
similar to the input IMF and for TV = 100 suffer from the 
low number of escaped stars. Velocity histograms are only 
shown for the cases with 100% initial binaries. All velocity 
histograms are also rather similar but the cases with bina- 
ries have a larger high-velocity tail when compared to the 
single star only run. 

4.1 Comparison with analytical estimates 

The Panels B of Figs. [3] and [7J show the number of stars 
lost per number of relaxations the clusters experiences dur- 
ing 5 Myr for TV = 1000 and TV = 100, respectively. Besides 
the results from the different initial conditions there are also 
shown the analytic estimates as given by eqs. [3] (solid black 
dots), [6] (open squares) and [8] (open triangles). In both cases 
the eq. [3] predicts far too low escape numbers which is ex- 
pected as this equation should be applied only to clusters 
with large TV. For the TV = 1000 case the resulting star loss is 
between the predictions of eqs. [S] and [5] but somewhat closer 
to the binary interaction dominated model (eq. [S}. This can 
be seen as an indication that star-loss in such clusters is not 
purely dominated by binary interactions (eq. but evapo- 
ration (eq.[8]) already plays a role. Though, the pure evapo- 
ration model of eq. |8]technically also only applies to clusters 
with larger TV. The sub-virial models are much higher above 
the predictions of eq. [6] probably because the shorter the 
time scales of the collapsing clusters lead to more binary 
interactions compared to the virial cases. 

The lower TV clusters with TV = 100 behave different in 
the sense that the numerical number of escapees is in be- 
tween the estimates given by eqs. [3] and [SJ but the binary 
interaction dominated model (eq. [6]) overestimates star-loss 
strongly, while the evaporation model (eq.|8]) seems to give a 
result closer to the numerical calculations. Here it might be 
that higher escape velocity due to the background poten- 
tial impedes the loss of stars after low-energy encounters. 
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As the radii for the N = 100 and the N = 1000 clusters 
are the same, the more massive clusters have higher den- 
sities and therefore higher encounter rates. For these it is 
therefore more likely that they experience high-energy en- 
counters compared to the N = 100 clusters and therefore the 
N = 1000 clusters might be less affected by the background 
potential. 

It should be noted here that none of the analytical esti- 
mates take any background potential into account and only 
focus on one physical process for star-loss, either energy 
equipartition or binary interactions. 
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APPENDIX A: MASS AND VELOCITY 
HISTOGRAMS 

In this appendix are shown the histograms of the mass and 
velocities of the escaped stars. Fig. I All shows a mass his- 
togram and velocity histogram example for iV = 1000 while 
Fig- EH shows the same histograms for N = 100. 



5 CONCLUSIONS 

• Low-mass star clusters loose between 1% and 20% of 
their stars (between 1% and 15% of their mass) within their 
first 5 Myr. 

• This is nearly independent whether or not the clusters 
have primordial binaries but collapsing clusters loose more 
stars. 

• The mean velocities of the escaping stars depend on 
the richness of the cluster and its radius and can be up to 
8.5 km/s for TV = 1000 and 7? ccl = 0.1 pc. This is slightly 
enlarged if primordial binaries are present. 

• The velocity dispersion of the escape velocities is about 
half of the mean escape velocity when less than 50% of the 
stars are in primordial binaries. In the case of high binary 
fractions it can be as large as the mean escape velocity itself. 

• The mean mass of the escapees depends on the richness 
of the cluster, its radius and the presence of binaries and 
vary up to an factor of 4. 

• None of the analytic estimates from the literature de- 
scribe the star loss precisely but they generally give the right 
order of magnitude of lost stars. 

• The currently known number of fast escapees in the 
ONC is consistent with the predictions of the TV-body cal- 
culations. 

Henceforth, as the majority of the escapees are not 
fast run-away stars but evaporate relatively slow, the pres- 
ence of 10% to 20% of distributed stars in a young star- 
forming region can be attributed to dynamical ejections 
from low- AT embedded star clusters. Especially, when tak- 
ing into account that un-resolved binaries and chance su- 
perpositions of stars will t end to underestimate the num- 
ber of stars within a c l uster (I Kouwenhoven fc de Griisll2008l ; 
iMafz Apellanizl 120081 . 120091 : I Weidner et alj|2009h . No such 
miscounting due to crowding is to be expected outside of 
the clusters, leading to a overestimation of the ratio of stars 
around clusters to stars inside of them. Current estimates 
for the amount of distributed stars in young star-forming 
regions are therefore not in contradiction with the assump- 
tions that most stars are formed in tight embe dded clusters , 
which is consistent with o ther JV-body studies l|Kroupal 19981 ; 
iKroupa fc Bouvierll2003T ). 
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Figure Al. The mass function {Panel A) of the lost stars for 
the models started without initial binaries and with a 1000 stars. 
Panel B: velocity histogram for the case with 100% initial bina- 
ries. All models started without initial binaries and 1000 stars. 
The solid histogram is the -R ec i = 0.1 pc model, the dotted his- 
togram is the R ec i = 0.2 pc model, the short-dashed histogram is 
the Reel = 0.3 pc model, the long-dashed histogram is the R e cl = 
0.4 pc model and the dashed-dotted histogram is the R ec i = 0.5 
pc model. The solid line in the mass function histogram is the 
input IMF (not to scale). 
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Figure A2. The mass function (Panel A) of the lost stars for 
the models started without initial binaries and with a 100 stars. 
Panel B: velocity histogram for the case with 100% initial bi- 
naries. All The solid histogram is the R cc \ = 0.1 pc model, the 
dotted histogram is the i? C cl = 0.25 pc model and the long-dashed 
histogram is the R ec i = 0.5 pc model. The solid line in the mass 
function histogram is the input IMF (not to scale). 
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